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We define the metamagnetic phase transition of itinerant electrons controlled by the spin injection 
mechanism. The current flow between a ferromagnetic metal and a metamagnetic metal produces 
the non-equilibrium shift of chemical potential for spin up and spin down electrons that acts as an 
effective magnetic field driving the metamagnetic transition. 

PACS numbers: 72.25.Hg, 75.30.Kz 



INTRODUCTION AND MAIN RESULTS 

The idea of the spin injection from ferromagneic metal 
to paramagnetic metal was first proposed by Aronov [l| . 
In the spin injection process the charge current flow be- 
tween the ferromagnetic and paramagnetic metals pro- 
duces the non-equilibrium magnetization in the param- 
agnet. This magnetization is proportional to the induced 
chemical potentials difference of electrons with opposite 
spins [lj - the spin accumulation. Non-equilibrium spin 
imbalance due to injection was observed by Johnson and 
Silsbee Q. The theory of spin injection was developed 
in details in many works [3|47| and well studied exper- 
imentally, see for a review [1, 0|. However, the degree 
of electron spin polarization is relatively small at stan- 
dard spin injection from ferromagnetic to paramagnetic 
metal jloL [Tl( ] . In order to increase the non-equilibrium 
polarization it is interesting to look for the possibility of 
spin-injection based magnetic transition in metamagnetic 
metals. Here we consider the metamagnetic transition of 
itinerant electrons induced by the spin injection mecha- 
nism. Let us briefly describe the properties of the meta- 
magnet [l2|, [l3| . When the energy splitting of electrons 
with opposite spins is smaller than the characteristic en- 
ergy scale of itinerant electrons, then magnetic part of the 
free energy density can be expanded in powers of mag- 
netization F(H, M) aM 2 + bM 4 + cM 4 - MH, where 
coefficients a, 6, c are determined by the energy depen- 
dence of the density of states at the Fermi level, H is the 
external magnetic field and M is the magnetization. 

At b < magnetic part of free energy F(H = 0, M) 
might have extremum at nonzero \M\ as it is shown in 
Fig. [TJ which schematically illustrates evolution of free 
energy with increasing magnetic field due to contribution 
of the term —MR. At small magnetic field the state 
with low magnetization has lower energy, while at mag- 
netic field larger than so called metamagnetic field H m 
the metamagnetic state acquires lower energy and system 
undergoes to state with higher magnetization. Metamag- 
netic state is induced by external magnetic field through 
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the first order phase transition 

Meta mag netic transition of itinerant electrons might 
appear 12, LL3| in strongly enhanced Pauli paramagnets 




when the Fermi level is close to peak in electron density of 



FIG. 1: Free energy F(H,M) dependence on the magnetiza- 
tion M of the metamagnet shown schematically for different 
magnetic fields H2 > H±. The state with higher magnetiza- 
tion has lower energy at higher magnetic fields. Inset: De- 
pendence of the magnetization on magnetic field. 



states. In this case Zeeman splitting increases the density 
of states and drives the ferromagnetic transition. 

The chemical potentials difference of electrons with op- 
posite spins is the analog of external magnetic field in the 
non-equilibrium case. Characteristic feature of this effec- 
tive magnetic field H e ^(x) is the spatial non- uniformity 
which results in the finite length of the metamagnetic 
state. In the ferromagnet - metamagnetic metal contact 
spin accumulation and therefore effective magnetic field 
is generated in the region of the order of spin relaxation 
length at the vicinity of contact with ferromagnet and 
at the domain wall between metamagnetic - paramag- 
netic states. We assume that the domain wall thick- 
ness is much smaller than the spin relaxation length. 
Schematically ferromagnetic metal - metamagnetic state 
contact is shown in Fig. [2j Metamagnetic state is lo- 
cated at < x < d. Metamagnetic phase emerges at 
electric currents such that the effective magnetic field is 
H eE (x = 0) > H m . If d is of order or larger than the 
spin relaxation length then effective field H eS (x) can be 
estimated as a sum H eE (x) = H^_ m (x) + H^_ p (x) of 
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the fields due to spin accumulation at boundary x = 
eJ 2R F R 



F-m 



(x) 
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and effective field due to spin accumulation at domain 
wall x = d 



gjlB Rm + ^Rp 



This case is shown by the solid line in Fig. 02 In expres- 
sions (pQ) and J is the current density, e is the electron 
charge, \ib = \e\h/2mc is the Bohr magneton and g = 2 
for electrons, 
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is proportional to the current polarizations, where <j a — 
e 2 D a v a are the corresponding conductivities in the fer- 
romagnetic, metamagnetic and paramagnetic states for 
electrons with spin a, D a is the diffusion coefficient, u a 
is the density of states at the Fermi level, 



RF,m — ^F,i 



n ^ = f (4) 



are the effective resistances and the spin relaxation 
lengths are defined as £ = \J Dt s , where in each state 
D = (D^a^ + D^ai)/(a^ + cr^) and t s is spin relaxation 
time. 

In case of small thickness of domain wall transition be- 
tween metamagnetic and paramagnetic states takes place 
at x = d when 



H eS (d) = H n 



(5) 



as shown in Fig. 03 Taking the sum of expressions (pp) and 
(|2|) we estimate the corresponding length of the metam- 
agnetic region at d > £ m as 



d. 
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At large electrical current when gjisHm/ eJ2R m — > 
RpU-m I (Rm + Rp) , according to expression (j6j) the length 
of metamagnetic region d — >• oo. Threshold current den- 
sity dependence of d occurs because of spin accumulation 
generation at domain wall between metamagnetic and 
paramagnetic states. However, since the effective field in 
most part of the region is below H mi the energy of sta- 
tionary metamagnetic state at large d becomes smaller 
than the energy of paramagnetic state. We propose that, 
system must undergo to paramagnetic state at large val- 
ues of current density. More detailed discussion of the 
transition is given below. 
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FIG. 2: Ferromagnetic (-L/2,0) metamagnetic (0,L/2) con- 
tact. Shaded region defines the high magnetization state of 
the metamagnet induced by the spin injection from the ferro- 
magnet. lF,£m,£ P are the spin diffusion lengths in ferromag- 
net and metamagnet in high and low magnetization regimes. 
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FIG. 3: Up: Dependence of the effective magnetic field on co- 
ordinate for two values of the current density | J2 1 > \Ji\. Ef- 
fective field decreases in the metamagnet and the phase tran- 
sition undergoes at x — d{J) when H eE = Hm- Down: Mag- 
netization profile, where Mi?, M m , M p are the corresponding 
magnetizations of the ferromagnet, high and low magnetiza- 
tion states of the metamagnet 



ELECTRICAL SPIN INJECTION 

Consider the spin injection process from the ferromag- 
netic metal to metamagnetic metal. We focus on the 
spin and charge transport in the presence of the spin- 
orbit coupling and the short-range exchange electron- 
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electron interactions. We assume the vector of the non- 
equilibrium magnetization in the metamagnetic metal to 
be parallel to the vector of the magnetization in the fer- 
romagnet. The Green's function in Keldysh technique in 
the matrix form 



G 



G R G K 
G A 



is given by retarded Cr (x, x'), advanced G A (x, x') and 
Keldysh function G K (x, x'), where x = (r,t) denote po- 
sition and time arguments, hat Q stands for the matrix 
in spin space. Further we will consider the stationary 
regime in which function G satisfies the equation 

[(e>+±VZ + H-U(r)-e<f>(r))i- 
-U so (r) + e(r)]G(r, v',oj) = 16(r - r') 

where 0(r) is the electrostatic potential, U (r) is the ran- 
dom potential of the impurities assumed to be Gaussian 
distributed with (U(r)) = 0, (U(t)U(t')) = 27ri/r<f(r-r'), 
r is the mean free time and v = {y^ + is the 

density of states. Spin-orbit interactions of electrons 
with impurities is described by the potential U so (r) = 
ija(\/U(r) x V), where 7 is the spin-orbit coupling con- 
stant and a is the Pauli matrix [14(. The contribution 
of the short-range electron-electron exchange interactions 
to the spin splitting is described by the term e(y) 



£a(r) 



-i\ f dpdu K 



where A is the electron coupling constant and it is con- 
venient to apply the Fourier transformation with respect 
to the relative coordinates 1*1 = r — r' as 



G(R,p,u;) = / ^riG(R + n/2,R-n/2)e" 



*pri 



in which R = (r + r , )/2. The retarded and advanced 
Green functions G R and G A averaged over disorder in 
the fir ^> 1 approximation are diagonal in spin and are 
given by 



■ £ P + e(R) ± 17 



(7) 



where £ p is the electron dispersion, 2j a 



tzi)/2 and r- 



te 1 



1 _ 
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r s J- is the inverse scattering 



times due to disorder and spin-orbit interactions for the 
spin state a, = 4/3r sa [1J]. We assume that the 
momentum relaxation time ro a is smaller than the time 
t sa corresponding to the spin flip process. Let us note, 
that we are considering the metamagnet when the ex- 
change energy is the coordinate dependent function. In 
deriving the equation for the function G K we assume the 
limit when the exchange energy is small compared to the 



Fermi energy | — \ / n < 1. In this limit the equation 
for the function G K yields 

v(V R + [V R e a + eE}d ep )G« = - - £ + j^) 



Fg—F- 

2t as 



here E = — V0(r) is the electric field and we denote the 
coordinate and frequency dependent function 

In the diffusion approximation for the function F a (R, uo) 
one obtains the equation 



VG a VF a (R,uj) 



F a (R,uj)-F_ a (-R,uj) 



2v 



(8) 



where G a = e 2 v a D a is the conductivity, D a = v^r a /3 is 
the diffusion coefficient and the density of states are the 
space dependent functions. 

Let us consider the system when functions in Eq. J5J) 
depend on one spatial coordinate (x) only. Consider the 
interface between a ferromagnetic metal that occupies 
the region (x < 0) and a metamagnetic metal (x > 0). 
We assume that the lengths of the ferromagnetic and 
metamagnetic regions L/2 are much larger than the cor- 
responding spin diffusion lengths. We also assume the 
external reservoirs of the sample at x = ±L/2 to be in 
the spin equilibrium state. The electric field in the sys- 
tem is treated through the boundary conditions 



F a (-L/2,u) = f(u - eV/2) 
F a (L/2,uj) = f(o J + eV/2) 



(9) 



where f(cj) = tanh(cj/2T) and V = E/L is the potential 
difference across the structure. The solution of the Eq. 
JHJ) is the continuous function at the interface x = 



F a (0-,uj) =F a (0+,cj) 
while the derivatives satisfy 



(10) 



dF a (x,u) _ dF a (x,uj) 

Cr* " \x=0- — — \x=0+ l iJ -J 



dx 



dx 



describing the continuity of the current density at the 
interface. The current density carried by spin up and 
spin down electrons is given as 



dF a (x,u) 
dx 



We solve Eq. (j8j) assuming the boundaries ferromagnet- 
paramagnet, ferromagnet-metamagnet and metamagnet- 
paramagnet independently. This approximation is valid 
in the limit when the length of the metamagnet d > £ m . 
Taking into account the length of the system to be larger 
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than the corresponding spin-diffusion lengths we solve 
the diffusion equation in the region x > with boundary 
([9]) and continuity ([TO]) , ([TT]) conditions. 



F tMPim {x, u) = f{u> + eV/2) + A p>m [(a; - L/2) ± 

P.m] (12) 



±2° l,t\p,m{HF - n p/m )— — ^ re 



(^F + Rp,m) 

where p, m denotes low and high magnetization regimes 
of the metamagnet and F stands for the ferromagnet, 
coefficient 

A 2(q tF + a iF )[f(uj + eV/2) - fju - eV/2)} 

is connected with the current density as 

J = Jf 0) + J±(x) = ^ J [ a t|p,m + Vl\p,rn]A p ,mdu 

The conductivity spin polarization and resistivities in the 
ferromagnet and metamagnet are defined by expressions 
dU and (|4j). Note, that in the low magnetization regime 
of the metamagnet <j^ p — cr^ p = <r p /2, D^ p = D^ p = D p . 
Solution ([T2]) has to be supplemented by the local neu- 
trality condition which self-consistently determines the 
electric potential in the sample. The spin injection pro- 
cess does not change concentration of electrons 



V = - / [v t F t (x, u) + v^x, u)]duj (13) 



PARAMAGNETIC STATE 

Low magnetization state of the metamagnet can be 
studied by solving Eq. J5J) assuming contact between 
ferromagnetic metal and paramagnetic metal at x = 0. 
Effective field due to spin accumulation is 



■eff/ \ e ^ 2RpR 



gfi B Rf + R P 
and magnetization at x > is 



1 — AVr> 



(14) 



Effective magnetic field produced by the spin accumula- 
tion in the ferromagnetic metal at x < is 



(15) 



here expressions for resistances Rf and R p are given by 
Eq. dU). 



METAMAGNETIC TRANSITION 

The self-consistency equation for the magnetization 
density M(x) in the sample is defined as 



(16) 



M(x) = g/2 B [si(x) - e t (x)}/\ = 
JWM x ^)-n F i(x,uj)]duj 



In the case of equilibrium metamagnetic metal, Eq. (jH 
has two solutions even without the external magnetic 
field, corresponding to two minima of free energy, see in- 
set in Fig. (pQ). Transition between these solutions takes 
place when magnetic field is equal to H m . One could 
verify that in linear on V response spin dependent part 
of expression ([T2]) enters Eqs. (p3J [16]) as magnetic field. 

The procedure of finding solutions is following. We as- 
sume that there is metamagnetic state in the system at 
< x < d. Then we solve Eq. (j8]) for the spin accumu- 
lation at two boundaries and self consistently determine 
the value of d from Eq. ([5]). 

To obtain Eq. ([5]) we need to consider transition in 
more detail. Near transition between metamagnetic and 
paramagnetic states we need to include the spatial deriva- 
tives of magnetization into consideration, so 

dx 2 5M 

Here K is positive constant. Let we have solution 
M w (x — d), describing transition between metamagnetic 
and paramagnetic states at point x = d in uniform mag- 
netic field. It is solution of equation 

d 2 SF{Hm ,M w ) 
- K d^ Mw + ^M~ w =° 

Assuming small difference H^(x) — H m at x ~ d and 
substituting M = M w (x — d) + JM, we obtain 

d 2 l S*F(H m ,M) s 
~ cbfi + 2 Jm 2 \m=mJM = H m (x)-H m 

Solution of this equation exists if 



dx^ (x)(H^(x) - H,, 







(17) 



where ^o(x) is eigenfunction, corresponding to zero Eq 
mode of equation 



V ^ T 

'dx 2 



1 5 2 F(H mi M) ] 
2 5M 2 



\M=M, 



E k 9 



k^k 



^oW has no zeros and is localized near x = d in region 
of order of domain wall thickness. This mode describes 
small translation of M, so in the uniform field one has 
Eo = 0. In the case when t m ^ p are much larger than 
domain wall thickness from condition ([TT]) one obtains 
Eq. ©. 
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Metamagnetic state 

Let us discuss different realizations of spin accumula- 
tion. 

1) . Consider the case when lip — H m and H m have the 
same sign. Effective fields of both contacts have same 
sign too. The estimation for the length of metamagnetic 
region d in the limit d > £ m is given by expression (j6j). d 
diverges at some threshold electrical current density. 

2) . Let Up — n m and U m have opposite signs. Thus, 
effective fields of both contacts have different signs too. 
Analysis shows that solution of Eq.(j5j) with finite d exist 
at \Hf_ m (0)\ > \H^_ p (d)\. With incresing electrical 
current density d stays finite. 

In metamagnetic region the magnetization is 

M m (x) =M° m + iit^Hfix) 

Here M m is the magnetization of metamagnetic 
state, calculated at zero magnetic field and v m = 

Zv^mVim/ (l/f m + Vim). 

Spin accumulation appears also in the ferromagnetic 
metal at x < as 

Hf m {x) = ^-^^[n F - n m ]e^ (18) 

gflB tip + K m 

here expressions for Rp and R m are given by (j4]). 



Free energy criterium 

We propose that the realization of metamagnetic state 
must be energetically favorable over realization of the 
paramagnetic state. In the linear on the applied current 
regime the magnetic part of free energy in the case of 
paramagnetic state realization is 



S^Fp = -M F 



f Hf p (x)dx 

J-L/2 



where Mp is the magnetization of ferromagnetic contact. 
In the case of metamagnetic state realization it is 



^Fmp - Mp J° L/2 *j. Fm 



H e JL(x)dx- 
H m ]dx + F s 



Effective magnetic fields in ferromagnetic region are given 
by expressions ([T5]) and ([15)1 . Fs is the energy, associated 
with domain wall and boundary F — m. While domain 
wall energy is positive, the sign of F — m boundary energy 
depends on relative directions of magnetizations in the 
ferromagnet and metamagnet. Estimation of Fs depends 
on details that are beyond the scope of the paper. 

From criterium JJ-pp — ^Fmp > for realization 
of the metamagnetic transition one can estimate the 



threshold value of current density. In the limiting case 
Rp > R m > R pi which assumes the contribution of the 
boundary m — p and F — p to the free energy is smaller 
than the corresponding contribution from the F — m in- 
terface, one can estimate as 



Jthr 



H rn /2R rj 



d(Jthr) 



e (n F -n m )(l+£ F M F /^ m M0) 



Ferromagnet-metamagnet-ferromagnet structure 

Let us briefly discuss the spin injected metamagnetic 
state in system with metamagnetic metal placed between 
two ferromagnetic contacts with opposite directions of 
magnetizations. In this case JJ-pp = 0, because of can- 
celation of contributions in ferromagnets with opposite 
magnetizations. In metamagnetic state 

«TmF = -M» foiH^(x) - H m ]dx 

Both ferromagnets contribute equally to the effective 
field. At d > £ m using expression (pQ) we obtain the 
threshold value of electrical current density at which 
(STFmF < as 



Jthr 



9^B 



Hrr 



Rp + R n 



d 



4RFR m [Up-U m } e n 



Note, that the expression for J t h r for the F — m contact 
in the limit discussed in the previous section is similar to 
the F — m — F contact. Also note, that the transition to 
paramagnetic state with increasing current is absent. 



CONCLUSIONS 

To conclude, we have studied the metamagnetic tran- 
sition of itinerant electrons in the metamagnet under the 
spin injection from the ferromagnetic metal. Spin in- 
jection produces the non-equilibrium effective magnetic 
field in metamagnet which drives the transition. We have 
calculated the effective magnetic fields and electrical cur- 
rents required for the metamagnetic transition. We have 
shown that the length of metamagnetic state has thresh- 
old dependence on electrical current due to the effective 
magnetic field self generated at domain wall. 

Typical values of the spin accumulation in metals are 
in the /ieV range [la, LL6[ , which corresponds to the ef- 
fective magnetic fields in tenth mT range at reasonably 
high current density. Metallic metamagnets with meta- 
magnetic field in tesla's range are well known [17]. Ap- 
plying external magnetic field one can easily bring such 
system close to the transition. 

Well studied YC02 in crystal form undergoes the meta- 
magnetic transition at H m = 70T [l3[, while in the 
nanoscale structured form it is a weak ferromagnet [181 ] , 
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suggesting the possibility of metamagnetic field reducing 
at proper technology. Other possibility is to study the 
system with temperature induced metamagnetic transi- 
tion [l9|. Unfortunately, spin relaxation length, the main 
parameter that governs the magnitude as well as the spa- 
tial distribution of effective field, is not known in such 
systems. 
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